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Wc obtain the double scaling limit of a set of superloop equations recently pro- 
posed to describe the coupling of two-dimensional supergravity to minimal super- 
conformal matter of type (2,4m). The continuum loop equations are described in 
terms of a c = 1 theory with a Z2- twisted scalar field and a Weyl-Majorana fermion 
in the Ramond sector. We have computed correlation functions in genus zero, one 
and partially in genus two. An integrable supersymmetric hierarchy describing our 
model has not yet been found. We present a heuristic argument showing that the 
purely bosonic part of our model is described by the KdV-hierarchy. 
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1. INTRODUCTION 



A description of the coupling of = 1 superconformal field theories [1] to 
two-dimensional supergravity in terms of discrete triangulations of super- Riemann 
surfaces is an interesting problem still unsolved. The continuum analysis carried 
out in [2] [3] has not yet found a complete discrete counterpart. The use of the 
Kadomtsev-Petviashvilii (KP) hierarchy initiated by Douglas [4] in the analysis of 
the double scaling hmit of the purely bosonic theories [5] was not fully succesfuUy 
pursued for the unitary superconformal (m, m -|- 2) models in [6] in terms of the 
Manin-Radul [7] supersymmetric extension of the KP-equations. The root of the 
problem was the incompatibility between the string equation and the fermionic 
fiows. In [8] a proposal was made to describe the coupling between (2,4m) su- 
perconformal models and world-sheet supergravity. The model is based on a set 
of superloop equations which are motivated by analogy with the bosonic case [9]. 
The planar solution to this model allowed the construction of bosonic and fermionic 
operators whose dressed gravitational dimensions are in agreement with the contin- 
uum hmit result for the (2, 4m)-models [2,3]. The planar solution was constructed 
for simplicity starting with even bosonic potentials. This as expected generated a 
doubling of degrees of freedom in the Neveu-Schwarz (NS) and Ramond (R) sec- 
tors of the theory. It was also possible to obtain the correlation functions of an 
arbitrary number of planar superloops. 

The guiding principle in [8] was a set of super- Virasoro constraints satisfied 
by the partition function, which code the superloop equations. Since the Virasoro 
constraints [10] [11] played a prominent role in the Witten-Kontsevich theory (the 
intersection theory of certain line bundles on the moduli space Aig,n of genus g- 
surfaces with n distinguished points) [12], it is reasonable to expect that a similar 
set of super- Virasoro constraints should capture some important features of the 
supermoduli space of super-Riemann surfaces. The model obtained in this way is 
an "eigenvalue" model. It is formulated in terms of a collection of even and A" 
odd eigenvalues (Aj, Oi). We have still not succeded in finding a description of the 
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model in terms of generalized matrices, whose large— expansion may provide a 
geometrical interpretation in terms of triangulated super-surfaces. The aim of this 
paper is to continue the analysis of this eigenvalue model. In section two we briefly 
formulate the model and its loop equations, in section three we solve it in the 
planar limit for general potentials. We identify the scaling limit, the critical points 
and see how the doubling of degrees of freedom is lifted for general potentials 
in section four. In section five we obtain the continuum limit for bosonic and 
fermionic loops. We prove, that in the continuum limit they become respectively 
a Z2-twisted scalar and a Weyl-Majorana fermion in the Ramond sector. Thus 
the super- Virasoro constraints in the double scaling limit are described by the 
super-energy-momentum tensor of a c = 1 superconformal field theory. In section 
six we obtain the solution to the continuum loop equations in genus zero, one 
and partially in genus two. For genus zero this provides a good verification of 
previous arguments. We find that the terms independent of fermionic couplings 
in the free energy coincides with the result of the bosonic model. This result is 
in shght conflict with the genus one calculation of [13]. However, there was no 
clear way of normalizing the supermoduli integration in [13], and therefore this 
may be the origin of the discrepancy. Motivated by the agreement up to genus two 
of the purely bosonic partition function of our model with the one appearing in 
the Kazakov multicritical points, we present a heuristic proof that this equivalence 
works to all orders of string pertubation theory. On the basis of our computations 
plus the heuristic proof, we conjecture that the KdV hierarchy describes the even 
flows of our model when all fermionic couplings are set to zero. This is indeed 
the hallmark for a supersymmetric extension of the KdV hierarchy, although at 
present we have not been able to identify it fully. In section seven we briefly review 
the difficulties we flnd in identifying our bosonic and fermionic flows with known 
supersymmetric extensions of the KdV- or KP-hierarchies [7,14,15,16]. In [8] the 
Mulase- Rabin [14] extension of KP seemed a promising candidate to describe our 
model. Recently some work in this direction has been done [17]*. The difficulties 

★ In this paper the double scahng hmit of the super- Virasoro constraints is also obtained. 
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we find may be related to the fact pointed out in [15] tliat tlie algebro-geometric 
solutions of the Mulase-Rabin or Manin-Radul hierarchies do not describe the 
moduli space of super-Riemann surfaces, but rather the moduli space of algebraic 
curves with genus g and a generic line bundle of degree g — 1- In this sense, 
our difficulties can be interpreted positively. It may be disappointing to find that 
the part of the free energy without fermion couplings coincides with the results 
of one-matrix models. Perhaps one should look at this result more positively, 
because it indicates that our model does describe geometric objects which extend 
fermionically Riemann surfaces, even though we have not formulated the model 
in terms of triangulations. The kinds of geometries described by the Super-loop 
Equations depend very much on the type of supersymmetric extension of the KdV 
hierarchy emboddied by our model. Section eight contains the conclusions and 
outlook. 



We first review some properties of the purely bosonic one-matrix models and 
see how they are generalized in the supersymmetric case. 

The general one-matrix model partition function is 



where $ is a Hermitian N x N matrix, and /i^ is the bare cosmological constant. 
The starting point of Kazakov's analysis of multicritical points [18] was a set of 
planar loop equations: 



2. SUPERLOOP EQUATIONS 






(2.1) 



A;>0 
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E ^9k-K7k^wil) = / dl'wil - r)w{l') , (2.2) 

k>i i 

w{l) describes a loop of lenght / bounding a surface with the topology of a 
disk. Equation (2.2) is the planar limit of the Schwinger- Dyson equations satisfied 
by the loop operator 



= = y (2.3) 

n=0 

in the general one- matrix model (2.1). Writing the partition function in terms 
of the free energy Z = ^, F = Fq + N~'^Fi + . . . the moments w^"") can be 
represented as: 

^(0) = A = -A2|^ , (2.4) 

dgn 

and the loop equations are equivalent to a set of Virasoro constraints satisfied by 
(2.1) [19]. They are obtained by implementing invariance of the partition function 
(2.1) under the change $ ^ $ + e$"^+^, n > -1 in (2.1): 

LnZ = n > -1 , (2.5) 



fy2 n q2 q 

^ 99n-kdgk ^ dgk+n 
Using the Laplace transform of the loop operator (2.3), 

oo (fc) 



e-Ml)dl = E ^ ' (2-7) 
and defining: 

pk+iqk+1 ' dgkdgi 



^(P^ ^) - E ;mTT Xk,i = A'^^ (2.8) 



k,l>0 

the loop equation, which is equivalent to the Virasoro constraints (2.5) and (2.6) 
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IS 

w{pf - V'{p)w{p) + -^x{P,P) = Polynomial(p) . (2.9) 

It is easy to check that (2.9) in the Umit iV^ oo becomes (2.2). The expUcit 
form of the polynomial on the right hand side is irrelevant. A useful way of rewrit- 
ing (2.9) is obtained after introducing an infinite set of creation and anihilation 
operators: 

a-n = --^ngn, n>0 ; = n>0, (2.10) 

AV2 N dgn 



and a scalar field: 



5<^(p) = E«np-"-' (2.11) 

with energy-momentum tensor: 

T{p) = ]^:d^{p)d^{p): (2.12) 

Then, (2.9) becomes 

Z-^T(p)Z = Polynomial(p) . (2.13) 

Furthermore, if we write the partition function (2.1) in terms of the eigenvalues of 
the matrix (Ai, A2, . . . , Ajv), but leaving the measure A^(A) undetermined : 

I nrfA,A2(A)e-xE.nA.) ^ (2.14) 
the constraints (2.5), (2.6) yield a differential equation satisfied by A 

whose solution up to a constant is the expected Van-der-Monde determinant 

^ = \\{\i-\j) . (2.16) 

In the supersymmetric case we proceed by analogy [8] with the above arguments 
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to obtain the corresponding loop equations. The loop operator depends on two 
variables (/, 9) where / is even and 9 is odd. We can define the Laplace transform 
as well 

oo 

w{p,n) = v{p) + Uu{p) = JdlJ d9e-P^-^^w{l,9) . (2.17) 



In particular, the Laplace transform of the operator D — ^ + is P = ^ — z-^. 
Asumming the loop w{l, 9) to behave well at Z = 0, oo, we can expand w{p, H) in 
inverse powers of p: 

-(^) = E^ ; -(^) = E^^' (2-18) 

A;>0 ^ fc>0 ^ 

v{p) and u{p) are respectively the fermionic and bosonic loops. To define the 
moments v^^\ u^^^ in terms of the free energy F = we introduce bosonic and 
fermionic oscillators: 



A (9 N 



a. 



a-p^--pgp , p= 0,1,2,... (2.19) 



P Ndgp ' A 



Ad . , , 

W/2 = ; b_p_y2 = -j-Cp+i/2 , p = 0,1,2,... (2.20) 

together with a free massless superfield: 

X{p,U) ^x{p) + Ui/j{p) (2.21) 

dX{p) = J2 «np-"-' m = E brp-'~'^^ , (2.22) 

n€Z reZ+1/2 
with the energy-momentum tensor: 

T{p, n) oc DXdX = ipdpx + U : {dpxdpX + dpipip) : (2.23) 
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The basic postulate in [8] is to take the superloop equations to be 

n)Z = Polynomial(p) . (2.24) 
In terms of u{p), v{p) these equations become: 

{u{p)-V\p)f + {v{p)-m\v{p)-ap)) + ^^^^^ + ^^^^^ , (2.25) 



where: 



u{p)v{p) - V'{p)v{p) - i{p)u{p) + = Qi , (2.26) 



k>0 k>0 



BF( \ ^ -^^ d'^F 



k,l>0 - 

^ i^oP'^'^i'^'^Skdgi (2-28) 

n>l r=l/2 ^ O-^r-C^n-r 

The quantities Qo, Qi are polynomials in p, and although their explicit form can be 
computed, they will not be needed throughout this paper. In terms of the original 
loop variables w{l,9) , the equations (2.25), (2.26)take a form similar to (2.2) 

VK.w(i,e) + 2K.Vw{i,e) = {woVw){i,e) , (2.29) 

with 

V = e-l^^ (2.30) 

d d^~^ 
/C = E(%fc^- 4-1/2)^ , (2.31) 

fc>i 

where the convolution between two superfunctions fi{z,9), f2{z,9) is defined ac- 
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cording to: 

z 

{fog){z,e) = J de' J f{z',e')g{z-z',d-d')dz' , (2.32) 


(see [8] for more details). The previous arguments suggest, in analogy with the 
one-matrix model, the introduction of a "superpotential" 

^) = E Y.^9k>H + efc+i/2^iA,') . (2.33) 

k>Q i=l 

The moments u^''^ and v^''^ (2-18) can thus be identified with derivatives of the free 
energy 

(^9n ?n+l/2 

Writing the partition function as: 

N 

Z = Y[dXid9iA{X,9)e-^^^^'^^ , (2.35) 

•' 1=1 

we can determine the explicit form of the measure A(A, 9) by imposing the super- 
Virasoro constraints (2.24). The explicit representation of the super- Virasoro op- 
erators using (2.19), (2.20) as differential operators is: 

(2.36) 

r=l/2 r=l/2 

Since {Gn-i/2) ^^7^-1/2} oc Ln+m-i, it suffices to implement — 0. This 
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leads to a set of equations: 



whose unique solution, up to a multiplicative constant is: 

A(A,^) = J](Ai-Aj-^i^,) . (2.39) 

Hence the model we would like to solve in the large- A?" hmit is: 

N 

H dXidOi Y[{Xi - Xj - ^i^,)e-x^(^.^) . (2.40) 

1=1 i<j 

The loop operator can be explicitly written as: 

w{l,9)^^Y.^^^^^''' ■ (2-41) 



Prom (2.40) and (2.41) one can derive (2.29). The simphfying assumption made 
in [8] was g2k+i = 0, A; > 0; i.e. the bosonic part of the potential (2.33) was 
taken to be even. In the next section we begin the analysis of (2.40) without this 
restriction. 
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3. SOLUTION TO THE PLANAR LOOP 
EQUATIONS: GENERAL POTENTIAL 

In this section we study the loop equation (2.25), (2.26) in genus zero for an 
arbitrary bosonic part of the super-potential. The simplifying assumption g2k+i = 
made in [8] generated a doubling of degrees of freedom in the Neveu-Schwarz and 
Ramond sectors of the theory that is not present in the continuum super-Liouville 
theory [20]. The planar loop equations follow from (2. 25), (2. 26): 

{u{p) - V'{p)f + (vip) - mYiviP) - e(p)) = Qo{p) , (3.1) 

(vip) - m)Hp) - V'ip)) = Qlip) . (3.2) 
Using the fact that Qi is fermionic, the solution to (3.1), (3.2) is: 

«W-V"(p) = v^-|^ (3.3) 

v{p)-aP) = ^0^ ■ (3.4) 

As in the pure gravity case we look for the one-cut solution. Since we make no 
assumptions concerning the parity of V{p), the one-cut solution takes the form: 

u{p) = uoip) + U2{p) = V'ip) - M{p)VA - ^ (3.5) 
v{p) = CiP) - ^ (3.6) 



with A = [p — x){p — y). The subindex in (3.5) indicates the order in fermionic 
couplings. We can also introduce variables R, S: 

y^S-yfR . (3.7) 
For V{p) — V{—p) the cut is symmetric and S — Since u{p) ~ 0(l/p), v{p) ~ 



10 



0{l/p) as \p\ oo, M{p),N{p) arc determined as functions of V'{p) and ^{p) 
respectively directly from this requirement. After M and are determined, the 
form of A{p) follows from demanding that the left hand sides of (3.1) and (3.2) 
must be polynomials. To write down the explicit form of M, N and A we note, 
that any analytic function /(p) can be written in the form: 



/(p) = /o(A) + A7i(A) 



A' ^^^2p-x-y^2(p-S) . (3.8) 
dp 

We split f(p) into two terms with opposite parity with respect to the change 
(p- S) ^ (S -p). Hence: 

M{p) = M"(A) + A'M+(A), 

N{p) = AN- (A) + A'iV+ ( A) , (3.9) 
A{p)^A-{A)+A'A+{A) . 

In N{p) we have used the fact that our solution should agree with the result 
obtained in [8], when V{p) ^V{-p) . The expansions of M, N and A in powers 
of A are given by: 



M±(A)=;^m±A^ A^± = ;^nfA^ A±(A)=5]a±A^ . (3.10) 

k>0 k>0 k>0 



To determine A, we substitute these expressions in (3.5), (3.6) and require the 
left-hand side of (3.1) and (3.2) to be polynomials in p. After some computations. 



-k To obtain the even case a factor of two must be taken into account as a normalization of 
the nps in (3.10) 
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we obtain that A (A) and A+(A) are completely given by Gq , Gq. The results are 

2R 



+ \2 



uIqUqUq — ARmQU^riQ) , (3-11) 



2R 

= — T^T—F^i'^o^i^Q-^Q^o^o) ■ (3-12) 



[m, 







4R{m+ 



Once we determine M{p) and N{p) we will have the complete solution to the planar 
model. Notice that N{p) will be linear in the fermionic couplings. This implies 
that the non- vanishing planar correlators contain at most two fermionic operators. 
This also occurs in the even case V{p) = V{—p), so that this rather surprising 
phenomenon does not depend on the type of potential chosen. We refer the reader 
to [8] for further discussions. 

The solution to the planar model given by eqs. (3.5) and (3.6) is parametrized 
by X and y. Since there is only one single physical parameter A, we should be able 
to express both x and y as functions of A. Let uo{p) be the purely bosonic part 
of the loop operator (no dependence on the couplings). Then the one-cut 

solution is: 

uo{p) = V'{p) - M{p)^{p-x){p-y) . (3.13) 
Since V'{p) and M{p) are polynomials in p, uq(p) admits the following expansion 

uoip) = C{x, y) + ^^^^ + 0(1) (3.14) 

p 

Dividing this expression by A-^/^ gives 

"0 V'{j)) 



M{p) 



_C{x,y) A(x,y) + ^C(x,y) 1 



P p^ P'' 



12 



which imphes 

Cix,y)^-{^- ^'^P^ (3.16) 
J 2m ^{p - x){p - y) 

oo 

A{x,y) ^ - f— (p — ) (3.17) 

oo 

But according to (2.18), uq{p) ~ 0(l/p). Thus we have the following constraint on 
X and y 

C(x,y) = (3.18) 

Equations (3.17) and (3.18) can be used, in principle, to rewrite the solution in 
terms of the single physical parameter A. These equations have appeared previ- 
ously in the matrix model literature (see for instance [21] and references therein) in 
a slightly different form. It is important to notice that the derivation is not based 
on the method of orthogonal polynomials [22] since this formalism is still missing 
in the super symmetric case. 

Note that the piece proportional to \{x -\- y) in (3.17) vanishes when (3.18) 
is enforced. Then one recovers the expression for A which is usually found in the 

literature. However, in order to compute partial derivatives it is essential to use 
the complete form (3.17). 

Expanding M in powers of p 

M{p) = Y,MnP^ (3.19) 



n>0 



we find 



and 



M„ ^ - ^ ^^p^ (3.20, 
'2m ^(q-x){q-y) 



J 2m(q-p)^{q-x){q-y) 
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In order to evaluate (3.11), (3.12), we also need explicit formulae for m^. These 
can be obtained as follows. Differenciating (3.17) and comparing the result with 
(3.21) gives 

d,A^^{x-y)M{x) , dyA^^{y-x)M{y) (3.22) 



On the other hand, by (3.10) 

M{x) = uiq + {x — y)mQ 
M{y) ^rriQ +{y- x)m^ 

This gives in terms of the single function A{x, y) 



^0 ^ (^^^)2 (^^ + ^l/)^ 



(3.23) 



(3.24) 



We conclude our analysis of the purely bosonic part of the solution by obtaining a 
simple expression for the derivative of uq with respect to the cosmological constant 
A. Under a variation {dx,dy) compatible with the constraint (3.16), duo is given 

by 



duoip) = -^[-^M{p)dA - MM{p)] 
dA 

p p^ 



(3.25) 



Since the expression multied by A is a polynomial, we must have 

dA = ^M{p)dA - AdM(p) (3.26) 



Therefore: 

duo{p) 



(3.27) 



dA V(p-x)(p-|/) 
When we discuss the continuum limit, the following loop operator will be rel- 
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evant 



uq{p) = uq{p) - V'{p) 



(3.28) 



Note that this also satisfies 

duo{p) 



(3.29) 



We now turn to the fermionic part of the solution. It is convenient to rewrite (3.6) 
in a slightly different form. Defining 



/c>0 



A' ~ 
v{p) + = i{p) - nVA (3.31) 



eq. (3.6) becomes 



where 

n^^-ll± ^(P) (3.32) 
° 2/ 27r^^(p-a;)(p-|/) 
oo 

and 

iV(p) = -/i|-— «IL== (3.33) 

oo 

We can also find an expression analogous to (3.27). Differentiating (3.31) with 
respect to x gives 

d,{v{p) + ^n+) = -^(^(p - y)N - Ad,N) (3.34) 

On the other hand, since by (2.18) v{p) ~ 0(l/p), 

d^ivip) + ^n+) = 2d^n+ + O(^) (3.35) 
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This implies 



l{p-y)N- Ada;N = a + 6A' (3.36) 



and if we set p = y we find a = {x — y)h. Since 



comparison with (3.35) yields h — dxTiQ . 
We finally obtain 



^ (a + 6A') = 26 + O(-) (3.37) 



/A 

Interchanging x and y gives 



^x(^^(p) + = -T (3.38) 



+ -7^ri+) = " ° (3.39) 



In order to complete the planar solution, we need explicit formulae for the functions 
Uq and appearing in (3. 11), (3. 12). First, note that eq. (3.36) can be written 

1 (p _ y)N - Ad^N = 2{p - y)dxn+ (3.40) 

Setting p — X and using (3.30) gives 

^{n^ + {x-y)n+)^2da:n+ (3.41) 



This, together with 



^inQ+{y-x)n+)^2dyn+ (3.42) 



(obtained from (3.39)) can be used to determine Uq and as functions of the 
single function Hq 

Uq = 2{dx + dy)n^ 

2 ^ , , (3.43) 



X y 



As in the bosonic case, we may define v{p) = v{p) — ^(p), which satisfies relations 
identical to (3.38) and (3.39). At this point we have determined the complete 
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solution to the planar model. In the next section we will consider the corresponding 
continuum limit. 



4. THE SCALING LIMIT 

The definition of the double scaling limit with an arbitrary potential involves 
some subtle considerations. These can be best understood by studying first the 
purely bosonic part of the theory. We take as the fundamental geometric 'observ- 
able' the macroscopic loop, defined by 



^ N 
u(l) = lim (V Af) = lim (4.1) 

i=l 



with n(P'l'^ — I fixed. Here m is a positive integer related to the order of criticahty. 
The continuum limit should be taken in such a way that (4.1) makes sense. This 
will be our guiding principle. Eq. (3.29) implies 

^ ^ 1 / dvv^ .^ 

^ > n-oo 27rz A 7 ^(p-a;)(p-|/) ^ ' 

c 

where the contour C encloses the cut (y, x). We may assume \y\ < \x\ < 1 without 
loss of generality. We will comment on the symmetric case y — —x later. Taking 
C as the unit circle around the origin with p — e^^ , 



27ri/ ^{p~ x){p-y) 27r/ ^(e^^ - a;)(e*^ - 7/) 



(4.3) 



We set Ac = 1, with A — 1 — aH, where t is the renormalized cosmological constant. 
X and y will approach the critical values Xc and yc, and N will be related to the 
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renormalized string coupling constant k by 

7V = a-2-V-(^,-^,)ll (4.4) 

(The factor {xc — Vc)^ is introduced for later convenience). Changing to the new 
variable z — a"^/"*^, eq. (4.2) becomes 



-dMl) = lim aV- v5^^± / ^"^"^ (4.5) 



a^O K 27r 



C 



This will vanish unless the integral itself is of order a-V"*, i.e., if X approaches 1 
as a^M, 

X = 1 - a^/'^u+ (4.6) 

Since x and y are not independent variables, y will approach its critical value yc 
(I yd < 1) a-t the same time 

y = yc + a'/"u- (4.7) 

Here n is another positive integer, which may be different from m. The integral 
in (4.5) is dominated by the region p ~ = 1, and the contour can be deformed 
into a straight line 

+ioo 

ndtUQ{l) = -^. I dz ^ ^ ' (4.8) 
2tti J y/z + u+ 

—ioo 

Here we recognize the definition of the inverse Laplace transform. Thus, if we 
define 

uq{z) = kjC[uo{1)] (4.9) 

we have 

dtuoi^) = -^2= (4.10) 

It is interesting to note that this result can also be obtained from (3.27) by a scaling 
p = 1 + a^/'^z in uo(p), together with (4.6) and (4.7). 
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There are two different ways of viewing uo{z): as the Laplace transform of the 
macroscopic loop uo{l), or as the continuum limit of the loop operator uo{p)- 

Prom (4.8) one gets the usual expression for the macroscopic loop 

dSoil) = rT (4.11) 

Note that uo{l) is independent of the scaling variable U-. This result holds as 
long as \yc\ < \xc\, independently of the values of m and n. It is easy to see that, 
for \yc\ = |.Tf :|. the dominant endpoint will be the one with the highest order of 
criticality. For a symmetric potential both endpoints contribute, and we get the 
familiar phenomenom of 'doubling'. As we shall see below, the situation is not so 
simple when one considers the fermionic contributions to the loop operators. 

In order to complete our description of the continuum limit for the purely 
bosonic part of the theory, we must consider the scaling of equations (3.16) and 
(3.17). Comparing first derivatives of A{x,y) and C{x,y), we find that they are 
not independent. Instead, 

^a;A ^l:{x- y)dxC 

I (4.12) 

Similarly, we have the following identities 

{x-y)dlyK^-\{d:,K-dyK) 

2 (4.13) 

{x - y)dlyC = -{d,C - dyC) 
Scahng (4.13) according to (4.6) and (4.7) gives 

dl_t ^ Aa^l''d+t + a^l'^d-t) (4.14) 

2(1 -yc) 

with an identical expression for C. For a — > the RHS vanishes, and we find 

d\_t = , d%_c = (4.15) 
where we have defined d± = and c = — |a~^(l — yc)C. Moreover, (4.12) 
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implies 

9+c = d+t , d-c = -d-t (4.16) 
The general solution to (4.15) and (4.16) can be written 

(4.17) 

c = --5^(t+<-V«-) = o 

where ip are the renormalized couplings. Eqs. (4.17) are the continuum version of 
(3.16) and (3.17). Adding and subtracting the equations in (4.17) give rise to two 
decoupled string equations 

^ = -5Z^>+' ^ = -$^V«^ (4.18) 

Comparing (4.17) and (4.18) we see that 

(f^t = (4.19) 



f± = --4^U±=o (4.20) 



1 

p\ 

The total derivatives in (4.20) are computed for variations du± consistent with the 
constraint (4.17). The relative factor of 2 between total and partial derivatives is 
important. It means that we can not set d-t = consistently. 

The m-multicritical point (at x) is defined by 

d+t\o ^ . . . ^ d'J^-Hlo ^ (4.21) 

One can impose similar constraints at y, with a different index n. But, as mentioned 
above, for \y\ < \x\ the continuum limit is controled by m (the converse is of course 
true for \y\ > \x\). 
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Eq. (4.18) can be used to compute the derivatives of u+ with respect to the 
renormahzed couphngs 

^ = uld,u^ , ^ = (4.22) 

These expressions will be useful in connection with the definition of the free energy. 

We now turn our attention to the fermionic contributions. The following iden- 
tity can be derived for Uq 

{x - y)d%n+ = ^id,n+ - dyn+) (4.23) 

and, if r is the scaling function corresponding to ^q", we have 

dl_T = (4.24) 

which implies 

^ = | + ^E(^;< + ^p-^"-) (4-25) 
p>Q 

This defines the renormahzed fermionic couplings f^. Using eqs. (3.5), (3.11), 
(3.12), (3.24), (3.43) we find the following expression for U2 

«2(p)~i(l-fc)^'V^-^""f^^ (4.26) 

2 (;2 + M_l_)'^/^ 



The fact that uq and U2 must scale in the same way fixes the scaling for 

1 - yc 



Then 



Mz) = 7^^^ (4.28) 



and we get for the bosonic loop 

^ ^ - ( 1 \m ) = - r^— 4.29 
In order to obtain the fermionic loop we first note that, in the continuum limit. 
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the RHS of (3.39) vanishes. This turns the partial derivative in (3.38) into a total 
derivative dx, and from (3.38) we can write 

ddv{z)+ - ^^ d^rd^ (4.30) 

where we have defined 

v{z)^a-^^l^"-^ (4.31) 

The different powers of a in the definitions of the continuum bosonic and fermionic 
loops are a consequence of their different scaling dimension. The macroscopic loops 
are obtained by the inverse Laplace transform. The result is 

2 IT , p-lu+ 
dtim - -X -e-'''^Td+Tdtu+) = ^ (4.32) 



1 , p-lu+ 

dt{v{l) + ^=Te-^"+) = 2^=d+Tdtu+ (4.33) 
nynl Kynl 

Note that the fermionic contributions to the loops arc not independent of 

In other words, even though the loops are dominated by the contribution from 

p ~ Xc, their expectation values depend on all the couplings, even those defining 

the behaviour of the potentials at y. This peculiarity of the supersymmetric theory 

implies that one has to be careful when trying to obtain the double scaling limit 

of the SuperVirasoro constraints. This will be our main concern in section 5. 

Multiloop correlators are obtained by acting on u{l) and v{l) with the appro- 
priate 'loop insertion operators'. We consider first the bosonic case. Since 

(,34) 

we have 

m-)--~^--)^j'm-) (4.35) 

We need to know the derivatives of x and y with respect to g^. These are obtained 
by differenciation of (3.16) and (3.17). Since we are interested in the limit A; — > oo. 
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we need the asymptotic behavior of integrals Uke (4.2). Converting the integral 
around the cut into an integral along the real axis, and using Laplace method for 
integrals dominated by endpoint contributions, we find 



X 

II 



dp p 1 X 



J{x-p){p-y) V^-y VTrk ' 
y 



\x\ > \y\ (4.36) 



and 



dx y^x-y fkx 2 dy i i ^ i i 

^ ~A — \/~^^' ^ 0' m>\y\ 4.37 

dgk 4 V TT da;A dgk 

These formulae can be used for an alternative derivation of the macroscopic loop. 
Scaling this expression gives the 'macroscopic loop insertion operator' 

It is important to note that this expression assumes that all the dependence on the 
couplings is given implicitly through u+ and U-. But one can always rewrite the 
expressions so that this is actually true. The Laplace transform of (4.38) is 

2{z + 'U+)3/2 Su+ 

This operator inserts a loop u{z) . 

The fermionic insertion operators are obtained along the same lines. The result 



IS 



and 



J^il) = -a-V- ^ (4 40) 



. 1 1 S 



J'{z)^-—^=— (4.41) 

It is easy to see that the macroscopic loops can be obtained by acting with 
and J^{1) on the following free energy 

n'^dlF = -u+ + 2dt{Td+Tdtu+ - rd-rdtu-) (4.42) 
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5. CONTINUUM SUPER- VIRASORO CONSTRAINTS 

Wc will now see that, in the continuum limit, the bosonic and fermionic loop 
operators become a Z2-twisted scalar bosonic field and a Weyl-Majorana fermion 
in the Ramond sector. This is the supersymmetric generalization of the well known 
result for the bosonic matrix model. However, proving this statement turns out 
to be rather subtle in our case. The reason is that, as mentionned in the last 
section, the correlators depend on all the couphngs, not just those describing the 
continuum limit at the dominant endpoint. It is not a priori clear that one can use 
a single Z2-twisted bosonic field and a single Weyl-Majorana fermion. However, 
this turns out to be true. 

First we will show that uq{z) can be identified with a free bosonic scalar field 
ip{z) in two dimensions with antiperiodic boundary conditions 

d^{z) = Yl «n+i/2^-"-'/' , ^{e'^^'z) = -<p{z) . (5.1) 
neZ 

The Laplace transform of the bosonic loop can be decomposed into two pieces 

u{z) = t{z) + u{z) (5.2) 

where 

k>0 ^ 

and t{z) is the non- universal part. t{z) admits an expansion in powers z'^~^^'^, for 
n > 0. 

Eq. (5.3) is in fact a definition of the scaling operators o"^. Since all the 
contributions to the loop come from p ~ = 1, it makes sense to expand in terms 
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of operators associated with the x-endpoint. Comparing (4.29) with (5.3) imphes 
d, (uiz) - f^±^l = - J__ + ^ (5.4) 
With the Laurent expansion of (5.4) for z — > oo: 

du{z) ^(-i)'=+ir(fc + |) /mJ+i k ^ ^ A 

= - > ^ T — -r^ h 2dt(uird+Tdtu+) z 

dt ^!r(i) \k + i + + t +)j 



-k-3/2 



the scahng operators (cr^) are given by: 



1 (-l)^r(A; + f) / M^+i 



(5.5) 



This result should also be obtained from the free energy (4.42) by differentiation 
with respect to f^. Then (5.5) imphes the foUowing relation between t]^ and t]^ 



t+ = ^ i . 5.6 

k ^ip^u k y ) 



We now turn to the computation of t[z). Since (2; + -U-i-)"^/^ can not contribute to 
the non-universal part, it is obvious that t{z) is a function of independent of 
the fermionic couplings. Then, subtracting (5.4) from (4.29) gives 



where we have used 



Thus t{z) satisfies 



t{z) = — ^ + (indep.of m+) (5.9) 
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On the other hand, according to (3.28), uo{z) has the following form 

uo{z) = -M{z)^/A = u'^^A''+'2 (5.10) 

n>0 

where M{z) is defined as the scaling hmit of M{p), and A — z-\- For — 0, 
uq{z) contains only positive half integral powers of z. This implies 

t{z)^UQ{z)\u+=Q--^ (5.11) 

Differentiating (5. 10) with respect to t and taking into account that dtUQ = — 1/ \/A, 
and using (5.7) yields: 

ul = 2d+t , d+u"^-^ = -{n + ^)«[J (5.12) 

and we find the following expression for ^0(2;) 

n>0 ^{'^+2) 

Setting = gives 

n>0 

and using (4.20), (5.6) we finally get 

t(z)=;^(A;+i)t+z^-V2 . (5.15) 

fc>0 

In total 

A;>0 A;>0 

The relations between the couphng constants of the model and the modes of the 
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bosonic field are 

d 1 t+ 

"n+l/2 = , "-n-1/2 = {^+7:)— n>0 (5.17) 

SO that we can write 

u{z) = KZ-^dip{z)Z (5.18) 

Wc proceed analogously with the fermionic loop. First, decompose into universal 
and non-universal parts 

v{z) = ri{z) + v{z) (5.19) 

with 



1 BF' 



n>0 ^ 



v{z) will be identified with a Weyl-Majorana fermion in the Ramond sector. The 
Laurent expansion of (4.30) determines the correlator of the scaling operators 



fc!r(|)_, 1 , A _ 2 , 



and the relation between and r/j is 



The computation of %]{£) follows closely the method used with t{z)^ and we 
simply quote the result 

"^^^ fc>0 

This makes it posssible to identify v{z) as a Weyl-Majorana fermion in the Ramond 
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sector. We can write 



v{z) = kZ-^iI){z)Z (5.24) 



In terms of the mode expansions (5.20) and (5.23), 



-n-1/2 



d T+ 



(5.25) 



while for the zero mode we have, 



2k otq 



guaranteeing V^q ~ I 

We have thus succeded in writing u{z) and v{z) entirely in terms of the cou- 
plings associated with a single endpoint. One can easily check that the loop inser- 
tion operators constructed in the last section can be written as 



J^(.) = ^.— (5.27) 



and 



J^(z)=5]z— ^-^ (5.28) 



This means that the couplings and t„ merely parametrize the loop correlators, 
but do not contribute to the 'dynamics'. 
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The super- Virasoro constraints in the continuum are therefore described by 
the super-energy momentum tensor of a single cT = 1 superconformal field theory, 

Tf{z) = ^dip{z)ij{z) , (5.29) 



Tb{z) = i : d^{z)d^{z) : +^ : d^P{z)^P{z) : -t-^ . (5.30) 



With the mode expansion 



^'h-)4e%5# . tb(^)=j:^2 . (5-31) 



n€Z neZ 



we obtain, in terms of the coupling constants: 



G 



^o^Or , Y^/K , 9 TO d ^ + d 

A;>0 n-hfe-hl " k>Q n+k+1 



^'^ fc>0 ^ ^Wfe ^ fe=l ^h-l^^n-h 
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They satisfy the conmutation relations 

{Gm, Gn} = 2Ln+m + \{^^ " ^) Wn,0 (5.34) 

In fact, the non-universal term 45^0 is fixed by these relations. 
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The basic postulate of [8] was to take the discrete superloop equations equiv- 
alent to: 



After the double scaling limit we obtain that the superloop equations in the con- 
tinuum are equivalent to: 



with T{z,9) = Tf{z) + eTB{z) given by (5.29) and (5.30) This proves one of the 
main results of this paper. The continuum limit of our superloop equations are 
described by a Z2-twisted massless scalar field, and a Weyl-Majorana fermion in 
the Ramond sector. 

We would hke to remark that the super Virasoro operators act on Z instead of 
\fZ, as happens for cases with doubling. This fact was conjectured in [10], but we 
do not know of any previous complete proof of this statement. 

We finish this section with a comment about the dimensions of scaling opera- 
tors. With a similar calculation as in [8] we obtain for (cr^), — k/m and for 
{v^) dk — k/m — l/2m. These are the gravitational scaling dimensions of opera- 
tors in the NS- resp. R- sector of (2, 4m) = 1 superconformal minimal models 
coupled to two dimensional gravity. Since we have considered generic potentials 
there appears no doubling of operators as it was in [8]. 



Z-^T{p, U)Z = Polynomial(p, U) 



(5.35) 



Z-'^T{z)Z = Polynomial(z) 



(5.36) 
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6. SUPERLOOP EQUATIONS IN THE CONTINUUM 



In the following we solve the continuum loop equations determined by (5.36). 
Here the formulae get more transparent than in the discrete theory. It is interesting 
to consider also the planar case and see how the results for the loop operators and 
the string equation found in section four appear in a simple way. For higher genera, 
i.e. for the torus, the double-torus, ... we obtain a sytematic expansion determining 
all correlators beyond the planar limit. In the arguments of this section we are 
going to take for simplicity t{u+,U-) as a function only of u+. In other words, 
we set — for all k. We first show how the loop equations are solved up to 
genus two, and then present a general heuristic argument showing that the purely 
bosonic part of our model is equivalent to the KdV-hierarchy. 

The two super loop equations in the double scaling limit, which are equivalent 
to the continuum super- Virasoro constraints are: 



i2{z)v{z) + K^x^^i^) = Polynomial (2;), 



(6.1) 



u{zf - v{z)dv{z) + K'^ix^^iz) + x^^{z) + A) = Polynomial(z), 



(6.2) 



where the two-loop operators are defined by: 




A;,/>0 



(6.3) 




31 



The loop operators and the free energy have the following genus expansion 

u{z) = Uo{z) + K^Ul{z) + . . . = ^ u'q'^\z) + K^uf^\z) + . . . , 

A;>0 

= v,{z) + k\,{z) + ... = J: vf^'\z) + kV^'^'\z) + ..., 

k>0 (6-4) 

X^'H^) = X^o\^) + «^'xS"^(^) + ■ ■ ■ a^BF,BB,FF , 
F = Fo + K^Fi + ... . 

The subindices in our notation indicate the genus, while for the order in fermionic 
couplings we introduce upper indices. 

Planar solution 

The leading terms in the genus expansion of equations (6.1) and (6.2) are 

uq{z)v(){z) — Polynomial (2;) , (6.5) 

uq{z)^ — vo(z)dzVo(z) — Po\ynomial(z) . (6.6) 
We follow closely the steps of the discrete case with a one-cut ansatz for u^^ (z) 



A{z) 



uo{z) = M{z)^/z + u + - 4^ , (6.7) 

[z + uY''^ 

vo{z) ^ . (6.8) 

Expanding in powers of {z + u) 

N{z)^Y,'^k{z + u)'' , 

k>0 

M{z) = ^mk{z + uf , (6.9) 

fe>0 

we see that A{z) is determined by demanding the right hand side of (6.5) to be 
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polynomial in z 



SoW^MWv/JTi;-^^-^. (6.10) 

M{z) is determined from w^^^ ~ 0(2;"^/^). We have: 

which holds only when 

lM{z) + {z + u)duM{z) = l^ . (6.12) 

Inserting (6.9) in this equation we obtain a relation between the coefficients mi^ 
and the renormalized cosmological constant Iq: 



ou 



rrik = 



9^0 

-If r(i) 



(6.13) 



2 r(/t + |) ^i^'^+i 

The modes of N{z) are determined by demanding 



vq{z) = - r]{z) = 0{z-'/^) . (6.14) 



This imphes ~ 0(-z~^/^), leading to 



^ = -ik+l)nk+i k>l . (6.15) 
Fermi-statistics and compatibility between the bosonic and fermionic loop gives 

m = -AduTiQ . 
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Thus 



2 noduHQ 



uo{z) = M{z)VJ+^ - — , , ";3" , (6.16) 

1 ^ 

^tt'ol^) = — ; — ■ — dtno . (6.17) 



With the ansatz 



and (6.14) we obtain 



Thus 



^0 = X] l^nU"' , 
n>0 



_ (-l)^r(n+i) 
" 2n!r(i) 



no = -r(M) , r(M) = ^f„M". (6.18) 
The final result coming from (6.16) and (6.17) is then 

TdtT 



uq{z) = M{z)^/z + u + 



(z + m)3/2 ' 



1 ^ 

^t^o(^) = ; — ■ — dT{u) , 
\/z + u 

and this coincides with the results of previous sections. 

It is nice to see that from the purely bosonic part of uo{z) we obtain the planar 
string equation, since from: 

= M{z) - = 0{z-') , (6.19) 

we obtain the expected result from the vanishing of the terms proportional to z~^: 

fe>0 ^ V2) 

Genus one solution 
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The two superloop equations obtained from (6.1) and (6.2) are : 

1 



2^01*1 - vidvQ - vodvi + Xo^ + Xo^ + ^ = Polynomial (2;) , (6.21) 



viUq + vqUi + Xo^ = Polynomial(2;) . (6.22) 

Wo and vq are already determined. The two-loop correlators are obtained from 
(6.3) for F = Fo 

BF 1 ( 1 ^ I ^ 

^° =-^(^T^l7^^*"^"' 
^° = 8 (,(IT^ - j - 2^^ (7T^ ' ^'-'^^ 



Xo 



1/1 1 



8 \{z + uY 



We solve now the equations (6.21) and (6.22) according to the order of fermionic 
couplings 

u, = Y,uf'^ v, = Y,vf'^'K (6.24) 

A;>0 A;>0 

For the fermionic loop operator we obtain 

^fc+1/2 



with 



DuDt 



★ In the following we will omit the z dependences 
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Where we have introduced the notation D = d/dto. The bosonic loop operator is 



fC — -L 

with the coefficients 

W3/2 = -pFT - 2^^^i/2 ^7/2 = 5rD(DuDr) 



12Du 

W5/2 = + 61^3/2^^ - 2tDD\ = rVy^Dr . 



(6.26) 



The value of V1/2 cannot be determined by the requirement that the left hand 
side of (6.22) is a polynomial in z. It follows from the consistency between vi{z) 
and ui{z). This condition is given by the requirement of commuting derivatives 

(6.27) 



dt^dTk dTkdto 
where Fi denotes the genus one contribution to the free energy. 

For the purely bosonic part of ui{z) (no presence of fermionic couplings) one 
easily sees that this is the same result as for the one-matrix model. The equation 
to be solved in this case is 

uo{z)ui{z) +Xo{z) + Polynomial(2;) , (6.28) 



1/1 1 

8 \{z + uy 



where uo{z) = M {z)y/z + u is the solution for genus zero. Equation (6.28) coincides 
with the order zero equation in fermionic couplings coming from (6.21) just because 
the bosonic two-point correlators (6.23) and (6.29) coincide in the bosonic part. 
Our result does not agree with the continuum Liouville calculation for genus one 
of [13]. The origin of this discrepancy may be the fact, that there was no simple 
way to normalize the supermoduli integration in [13]. 
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The reader may have noticed, that the solution presented for vi{z) is only 
linear in the fermionic couplings, as it was for genus zero. This comes from the 
fact that, the third order equation in fermionic couplings 

v^^u'q^ + wqW]^^'' + v\^^ ^3/2 ~ Polynomial(2;) (6.30) 

(3) 

is solved hy v\ =0 which means that also on the torus we have a maximal 
coupling of two fermions. 

Genus two solution 

We consider now the situation in genus two to see if the free energy changes 
with respect to the one-matrix model. The two superloop equations to be solved 
are 

2uqU2 + Ui — V2dvQ — vodv2 — vidvi + Xi^ + Xi^ — Polynomial(2;) , (6.31) 



voU2 + uoV2 + uivi + Xi^ = Polynomial(2;) . (6.32) 

To calculate the free energy, we are only interested in the purely bosonic piece of 
the two loop correlators, which are given by (take (6.3) for F — Fi) 



FF ^ bb(0) ^ 13 1 3 mi 1 /3m^ 5 m2\ 1 . . 



We see, that again these two loop correlators coincide with the one-matrix model 
values. This immediately implies, that also in genus two we have the same value 
in the purely bosonic part of the free energy as in the one-matrix model. The 
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complete solution to order zero in fermionic couplings is 



"2 ^3/2 ^ ^5/2 ^ ^7/2 ^ ^9/2 ^ ^11/2 

GSm^ l^m\mi 145m2 llm\mz ^ 105m4 



32m^ 16m^ 128m^ 32m^ 128m^ 

^ GSmf 87TOim2 105m3 
~ 32m^ 32m^ ^ 128m^ 

63m? 145m2 , 203mi 105 



(6.34) 



32m^ 128m^ 128m^ 128mg 

We now continue the calculation for the case of pure gravity (mfc = 0, k >2). Our 
results reproduce, up to genus two the values expected for the Painleve-I equation. 
Since 

1 2 

mo = 7— mi = -rrinD^u (6.35) 
Du 3 

we get for the second derivative of the free energy the expansion 

M = - 4 + (^J - j + . . . . (6.36) 

For pure gravity we take u — y/t, to obtain 

which agrees with the first three terms of the solution to the Painleve-I equation 
appearing in pure gravity [5]. 

General properties of the Free Energy 

We now give the heuristic argument showing that the bosonic loop operator 
u{z) coincides with the pure gravity case with an even potential when all the 
fermion couplings are set to zero. First we write the loop equations for the pure 
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gravity case derived from an even potential. They are given by 

^u'^ + K'^ixiz) + ^) = Polynomial(^). (6.38) 

Note that equation (6.29) imphes xi^) + ~ where A — z + u. The 

correlator xi^) can be written in terms of the loop creation operator as in the 
supersymmetric case 

x{z)^Jf{z)u{z). (6.39) 

The explicit form of the loop insertion operators can be found in (4. 41), (5. 28). 
Since uq — M{z)A2^ Un{z) for n > will be obtained in terms of the negative 
half-integer powers of A, 

Un{z) = Y,<^~''~'^^- (6-40) 

k>l 

The expansion in (6.40) starts with A~^/^ because of the general structure of the 
loop operator discussed previously. Similarly, by (4.34) Jb ~ 0(A~^/^), and the 
first power of A in x{z) is A~^. These general remarks about the power series 
representation of x(2;0,m(2;) apply beyond genus zero. On the sphere there are 
non-universal contributions like 1/z'^ in (6.29) and we do not have such simple 
expansions in A. These considerations will be important in the following. 

The supersymmetric loop equations are given by (6.1,2). Since we are only in- 
terested in the leading contribution u^^\z) i.e. to order zero in fermionic couplings, 
it suffices to keep x^^ ^iid x^^ to order zero and v{z) to first order in fermionic 
couplings. Similarly, we can drop the term vdv in (6.2) which becomes, 

u{zf + K^{x^^{z) + x^^(^) + ^) = Polynomial(z). (6.41) 

We showed by direct computation that Xn^{z) — Xn^{z) at least for n = 0, 1, the 
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sphere and the torus. If this were true to all orders, then (6.41) would become: 

^u{zf + k\x^^{z) + ^) = Polynomial(^), (6.42) 

which is identical to (6.38), and our proposition would be proved. We will argue 
by induction that indeed Xn^i^) = Xn^i^)- Assume this to be true up to order 
n — 1. Expand (6.41) in powers of k?. This yields, 

n-l 

2uoUn + J2 '^kUn-k + = ^'(1), (6.43) 

A;=l 

where we have used Xn-ii^) ~ Xn-i(^)' have labeled a polynomial in z as 

C(l) in A. We have included the non-universal factor t{z) in uq. Expanding (6.2) 
in powers of we obtain 

n-l 

UQVn + UnVo + ^ UkVn-k + Xn-1 = ^(1) ■ (6-44) 
A;=l 

Acting with Jb on (6.43) yields 

n-l 

2uoxr + 2xr«n + 2 J2 ^n-kxS-k + 2 JsXn-i = 0(A-3/2), (6.45) 

k=l 

next, acting with d^Jp on (6.44) we obtain 

n-l 

Uoxn^) + «nXr + + dzJFXn-1 = ^(A"'/'), (6.46) 

k=l 

where we have used the induction hypothesis. For the same reason, we have the 
equalities 

d^JpXn-l = d,JJBVn-l = JsXn-l = ^SXn-1, (6-47) 

and we find that (6.45) and (6.46) are identical. Since uq — 0(A~^/^), this would 
imply that Xn^ — Xn^: ^ was to be shown. The basic subtlety in this argument 
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has to do with the equahty between Xq ^'Hd Xq ■ This is because naively these 
quantities contain infinities. To regulate them we have to point-split the variable 
in the loop creation operator with respect to the variable z in which we express the 
loop equations. This makes the count of powers of A rather subtle, and to assure 
that no mistakes are made we have to study in detail the two-loop equations. 
Nevertheless, we beheve that after the appropriate subtleties of the mathematical 
details are clarified, the conclusion of this heuristic proof will remain true. 

7. RELATION TO SUPERINTEGRABLE HIERARCHIES 

Up to this point we have been able to calculate correlation functions for the 
first cases in the genus expansion. Although this method is straightforward to 
apply for any genus, the procedure is cumbersome and of course it is essential to 
find the relation to supersymmetric integrable hierarchies. 

This relation is well established for the one-matrix model [4, 5, 23] . If we 
denote the specific heat with U = (PP) where P is the puncture operator, then 
the different multicritical points are characterized by the string equation 

oo 

-X = J2tnRn[U], (7.1) 

n=l 

here Rn[U] are the Gelfand-Dikii polynomials of the KdV-hierarchy, defined through| 
the recursion relations 

DRn+i[U] = {k^D^ + AUD + 2{DU))Rn[U] (7.2) 

The renormalization group flows of the models described by (7.1) coincide with 
the flows of the KdV-hierarchy: 

^^^DRn+^[U] (7.3) 

(7.2) implies therefore recursion relations between the flows. Equations (7.1) and 
(7.2) contain all the information about the correlation functions of the model. 
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Through the recent work of Witten and Kontsevich [12] we can directly relate 
the Virasoro constraints L„r = 0, n > — 1 with the r— function for the KdV 
hierarchy with initial condition L_ir = 0. Since in the supercase we do not yet 
have a formulation in terms of generalized matrices, it is useful to take a more 
pedestrian approach and see how the KdV-hierarchy emerges from the explicit 
solution of the loop equations. In the purely bosonic model, the one-point functions 
are given by: 

The definition of U as the two-point function of the puncture operator implies: 

D{ao)o = C/(°) = -\u (7.5) 
Using the planar string equation and (7.5) we learn that 



At genus one. 



^(i)r^) = ^ I I I (7 7) 

^ ^ 8ml{z + uf/^ 8mo(^ + «)5/2 ^ " ^ 



2 (-l)"r(n + f) f mi . n 
{an)i = K — 7^ — u + -u (7.8) 



Therefore, 

2mon!r(i) \2mo 3 
Hence the "heat capacity" U is given to genus one by: 

Similarly, to this order we can compute dU/dti, and obtain the well known KdV 
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equation 



— = K^D^U + 3DU^ 



(7.10) 



Hence, to this order: 



R2[U] ^ K^D^U + 3U^, 



(7.11) 



and 



Rn+l[U] 



22"+ir(n+ I) 




n!r(i) \n + l 



Including now the genus two correction to U which was computed in the previous 
section: 



it is easy to see that at genus two there is no correction to (7.10), and furthermore, 
that to this order wc have agreement with (7.2). Ahhough this is by no means 
the cleanest way to exhibit the equivalence between KdV and the loop equations, 
it is clear that to determine the differential operator appearing on the right hand 
side of (7.2) (assuming its existence) we need to know explicit correlators only for 
genera zero, one and two. Equation (7.10) together with the fact that it is not 
corrected in genus two and the perturbative understanding of (7.2) provide strong 
hints that the flows in the one-matrix model are governed by the KdV-hierarchy. 
In the supesymmetric case, we could proceed by analogy with the previous argu- 
ments. Since we have computed all correlation functions for genus zero and one, 
and we have partial results in genus two, we can try to explore what type of dif- 
ferential relations allow us to express the flows d/dtn, d/drn in terms of the basic 
flows d/dto, d/dTQ. We can also examine the supersymmetric extensions of the 
KdV- or KP-hierarchies. For KdV, the papers [16] find a one parameter family of 



U = C/(0) + «2[;(1) + ^4^(2) 



(7.13) 
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supersymmetric extensions given by: 



u = -u'" + 6uu' + Sia^f, 



(7.14) 



aC'" + (2 + a)uC' + 3u'C, 



(7.15) 




u (rcsp. ^) is a bosonic (rcsp. fermionic) function. This system is integrable only 
for a = 1,4. Only the case a = 1 is really supersymmetric. If we formally set 
^ = 0, we obtain KdV. Equations (7.14) and (7.15) can be therefore considered 
as supersymmetric extensions to the KdV equation. If on the other hand, we are 
also interested in fermionic flows {d/dTn), we need to extend these systems to 
incorporate odd flows. This has been done in [7, 14, 15]. In reference [7] the KP- 
hierarchy is extended using a supersymmetric Lax-pair formulation. A different 
approach is taken in [14, 15], where only the even (i.e. bosonic) flows admit a 
Lax pair representation. The fermionic flows have a simpler expression than those 
of [7]. A common feature of all these supersymmetric hierarchies is that if we 
consider only even flows, and we set all the fermionic variables to zero, we recover 
the standard KdV- or KP-hicrarchies. In terms of loop equations this is similar to 
formally considering only the bosonic part of the energy momentum tensor Tb{z) 
and ignore the fermion held dependence. This would generate the KdV-hierarchy 
as in the bosonic case. In this respect our loop equations are a good starting point 
to construct another supersymmetric extension of KdV. 

If we take the fermionic dependence into account, we have not yet been able 
to identify the corresponding differential relations which should lead to a local 
integrable hierarchy. Our analysis does not seem to make our one- or higher-point 
functions compatible with known extensions (7.14) and (7.15) or those in [16, 14, 
15]. This is perhaps not so bad as it sounds. If our model indeed describes a non- 
critical superstring, the alleged hierarchy describing the continuum limit should 
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capture the geometry of the supermoduU space of super-Riemann surfaces. From 
the work in [15] we know that the algebro-geometric data needed to construct quasi- 
periodic solutions to the Manin-Radul or Mulase- Rabin hierarchies does not include 
super-Riemann surfaces. We continue looking for an integrable supersymmetric 
hierarchy compatible with our loop equations. 

8. CONCLUSIONS AND OUTLOOK 

In this paper we have taken the double scaling limit of the superloop equations 
proposed in [8]. Working with general potentials, we have shown that the spectrum 
of anomalous dimensions coincides with those which follow from [2,3], when one 
couples two-dimensional supergravity to minimal superconformal matter of type 
(2,4m). The continuum limit of these theories is described by a Z2-twisted scalar 
representing bosonic loops, and a Weyl-Majorana fermion in the Ramond sector 
representing fermionic loops. We have solved the continuum superloop equations 
in genus zero and one completely, and partially in genus two. The piece of the free 
energy independent of the fermionic couplings is given in the pure supergravity 
case by the same expression as for the one-matrix models [5] i.e. by the same 
solution to the Painleve-I equation. This conclusion was shown to hold in general 
by a heuristic argument in Section seven. So far we have not been able to iden- 
tify an integrable superhierarchy, which reproduces our correlation functions. As 
explained in the text, this is not necessarily negative. No superhierarchy is yet 
known, which incorporates fully the geometry of super-Riemann surfaces. Finally, 
the corresponding generalization of multimatrix models in our context is still miss- 
ing. Work in these directions is in progress, and we hope to report on the results 
elsewhere. 
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